Mirages, ant i- mirages, and further surprises in quantum corrals with non-magnetic 

impurities 
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We investigate the local density of states (LDOS) for non-interacting electrons in a hard-wall ellipse 
in the presence of a single non-magnetic scattering center. Using a T-matrix analysis we calculate 
the local Green's function and observe a variety of quantum mirage effects for different impurity 
positions. Locating the impurity near positions with LDOS maxima for the impurity free corral can 
either lead to a reduction or an ehancement of the LDOS at the mirror image point, i.e. a mirage 
or anti-mirage effect, or even suppress LDOS maxima in the entire area of the corral. 
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I. INTRODUCTION 

The technological advances in scanning tunneling mi- 
croscopy (STM) have made it possible to manipulate in- 
dividual atoms on metallic surfaces [1]. These remark- 
able achievements allow local measurements of electronic 
properties with a spatial resolution of atomic length 
scales even on artificially designed geometries for surface 
adatoms [2]. The experimental advances have led e.g. to 
the spectacular observation of mirage effects in elliptic 
quantum corrals of magnetic Co atoms on Cu (111) sur- 
faces [3] . If an additional Co atom is placed at one of the 
foci of the ellipse, a Kondo resonance in the local density 
of states (LDOS) is not only observed at the magnetic 
ion itself but at the other, impurity free focus as well. 
The mirage effect can be viewed as a beautiful manifes- 
tation of quantum mechanical interference phenomena as 
a result of the multiple scattering events of the electrons 
from the impurity and the atoms forming the boundary 
of the corral and from the impurity at one focus. 

The theoretical work on this problem emphasised the 
many-body Kondo physics of itinerant electrons in the 
corral interacting with the localized magnetic moments 
of the impurities [4-7] . But it was also found that for the 
specific elliptic geometry even stronger mirage effects can 
be observed, if the additional impurity is moved slightly 
away from the focal point [6]. Furthermore, mirage ef- 
fects may even be absent depending on the intrinsic level 
width [5] or for specific electron densities inside the el- 
liptic corrals; one of the crucial quantities is the LDOS 
at the position where the additional adatom is placed. 
Therefore, already without invoking the Kondo physics 
the unique geometry of an elliptic corral appears to give 
rise to intriguing quantum mechanical effects which de- 
serve a detailed analysis also for the significantly simpler 
problem of a non-magnetic impurity in a hard-wall corral. 

Indeed we will show in this paper that a surprising 
variety of structures in the LDOS can be generated de- 



pending on the precise location of a single additional non- 
magnetic scattering center inside a hard-wall elliptic cor- 
ral. At the mirror position of the local impurity potential 
the LDOS may be suppressed or even enhanced, and also 
the almost complete suppression of LDOS maxima can 
be achieved. These observations underline the richness 
of quantum mechanical interference phenomena in an el- 
liptical geometry which may be verified experimentally. 



II. SOLUTION OF THE 2D SCHRODINGER 
EQUATION IN A HARD WALL ELLIPSE 

In a first step we review the solution of the 
two-dimensional (2D) Schrodinger equation for non- 
interacting electrons in a hard wall ellipse. Following 
Ref. [8] it is convenient to introduce elliptical coordinates 
through the transformation [9] 



x = aecos(9) cosh(ry) , 
y = ae sin(0) sinh(?7) , 



(1) 



where a denotes the semimajor axis and e is the eccentric- 
ity of the ellipse. In the new coordinates the Schrodinger 
equation takes the following form: 



2m(ae) 2 cosh(2?7) - cos(26<) 
V(r l )U(0,r,)=e^(e,r 1 ) 



de 2 



drj 2 



(2) 



where m is the electron mass and e is the energy eigen- 
value. Due to the hard-wall condition the potential V 
vanishes inside the elliptic corral and is infinite other- 
wise. Using the factorized ansatz for the eigenfunctions 



(3) 



the Schrodinger equation inside the ellipse (V(tj) = 0) is 
rearranged as 
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2kcos(29) 



0(6) 
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dr] 



2 + 2/ccosh(27y) 



$(»7) 
(4) 



where fe = (ae) 2 me/2h 2 . The Schrodinger equation thus 
separates with respect to the elliptic coordinates 9 and r\ 
and reduces to the two differential equations 



do 2 
d 2 

drj 



0(9) + (a-2k cos(26>)) 9(0) = , 



(5) 



^2 $(7?) - (a -2k cosh(2r?)) $(rj) = (6) 



where a is a separation constant. Equations (5) and 
(6) are the Mathieu equation and the modified Mathieu 
equation, respectively [9,10]. 

Since 9 is the polar angle and thus 0(9 + 2it) = 0(9) 
we select the periodic solutions of Eq. (5), which are the 
Mathieu functions of the first kind of integral order [9] . 
We obtain two types of solutions: 



or 



O r (9) = ce r (9,k c 
O r (6)=se r (6,k s 



(7) 
(8) 



where "ce" and "se" are the abbreviations for "cosine- 
elliptic" and "sine-elliptic" introduced by Whittaker [9]. 
Due to the periodicity condition there exists only a dis- 
crete set of numbers a — a r (k) for fixed k [10]. The index 
r in Eqs. (7) and (8) denotes the order of the Mathieu 
functions of the first kind. 

The solutions of Eq. (6), which are necessarily re- 
stricted to the same set of numbers a r (k) as the cor- 
responding solution of (5), are 



or 



M»7) 
$ r (r?) 



Ce r (r), k c 
Se r (rj, k s 



(9) 
(10) 



Ce and Se in Eqs. (9) and (10) denote the modified 
Mathieu functions of the first kind of integral order [9] . 
They must meet the hard- wall condition $>r(vo) = at 
the boundary line of the ellipse, which leads to a discrete 
set of values k n . This implies that the discrete sets for 
kZ and kt must meet the conditions 



or 



M%) = Ce r (r?o,££) , 
3v(r?o) = Se r (ij ,K) . 



(11) 
(12) 



Therefore, k^ and fc* are the n th zeroes of Ce r (r]o, k c ) 
and Se r (r]o,k s ), respectively, for fixed r. rjo is directly 
related to the eccentricity e of the ellipse through 



e = 



cosh(ryo) 



(13) 



Altogether we find the following form of the exact eigen- 
states inside the hard- wall ellipse [8]: 



VV C ,„ C (M) = ce r (9,k c n )Ce r (rj,k c n ) , 
V>r,n s (0, v) = se r (6, k s n )Se r (r), k s n ) . 



(14) 
(15) 



r, n and c(s) enumerate the quantum numbers for the 
eigenstates; the eigenenergies are determined by the 
hard-wall condition and the periodicity requirement for 
the Mathieu functions. Note that ce, se and Ce, Se are 
all real. The symmetry under reflection at the semimajor 
axis is even for the ce r function (cosine-elliptic) and odd 
for the se r function (sine-elliptic) for each r. Moreover 
the parity of ce r and se r is even/odd, if r is even/odd. 
Specifically, 



#(-M) = #(M) . 

r r (-e,v) = -r r (o,v) • 

r 2 i s (e + n, v ) = i>t(e, n ), 



V> 2 %(M). 



(16) 
(17) 
(18) 
(19) 



In the following we will specifically investigate an el- 
lipse with eccentricity e = 0.5. The energy scale for the 
eigenenergies is furthermore determined by the size of 
the ellipse, i.e. the length of the semimajor axis a. For a 
specific choice for a we orient ourselves at the typical size 
of the elliptic corral of Co atoms on a Cu (111) surface 
used in the experiments of Manoharan et al. in Ref. [3] 
for which a ~ 71 A. In this setup the Fermi energy tp 
is 450 meV which corresponds to the energy of the A2 nd 
eigenstate, and due to the spin degeneracy amounts to a 
particle number of 84 electrons inside the ellipse. For the 
specific choice for the size of the ellipse the ground state 
energy 



E = k 



2h 2 



(ae) 2 m 



(20) 



is Eq = 13 meV. For fixed eccentricity e = 0.5 and 
k n — &42 the wavefunction of the eigenstate at the Fermi 
surface has a high probability density at the two foci of 
the ellipse - a precondiction to observe strong mirage ef- 
fects when an additional impurity atom is placed at one 
of the focal points. 



III. THE SCATTERING PROBLEM 

The LDOS of the non-interacting electron system in- 
side the corral is easily obtained from the retarded one- 
particle Green's function 



GS et (r,r', £ ) = E 



V-;(r)V*(r') 



€j + : 



(21) 



where the eigenfunctions tpj are given by Eqs. (14) and 
(15) and ej is the corresponding eigenenergy of the ) th 
eigenstate. The free LDOS N (r, e) then follows from 
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FIG. 1. LDOS No(t, €42), £42 = e_p, for non-interacting 
electrons inside an ellipse with semimajor axis a (horizon- 
tal) and semiminor axis b (vertical). This viewpoint is the 
same for all plots in Figs. 1-4. The scale of the LDOS can 
be taken from Figs. 5-7. 
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FIG. 2. LDOS JV(r, £42) with an impurity (•) at the 
right focus ro = (ea, 0); U = 16.6-Eo^l (A = nab is the area 
of the ellipse). 
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FIG. 3. LDOS 7V(r, £42) with an impurity (•) at FIG. 4. LDOS N(r, £42) with an impurity (•) at 

r = (i, 0) = (0.22a, 0); U = 16.6E A. r = (i, 0) = (0, 0.33a); U = 16.6E A. 
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N (r,e) 



— ImGo(r, r,e) . 

7T 



(22) 



In Fig. 1 we present the LDOS at the Fermi energy 
ep = e 42 using a finite broadening 5 = (Ae) = 0.77 E 
where (Ae) is the average level spacing of the eigenenergy 
spectrum (S = 10 mcV for the parameters mentioned 
above for the ellipse with a = 71 A). While the eigenstates 
of the hard- wall ellipse are naturally sharp (5 = 0), in 
actual experiments the lifetime of the surface electrons 
is finite e.g. because of the imperfect hard-wall condi- 
tion. Using a finite broadening is therefore reasonable 
not only for practical purposes. We note that iV (r, e 42 ) 
is dominated by the contributions of the 42 nd and 43 rd 
state, because the energy difference between these states 
is seven times smaller than the energy difference between 
41 s * and 42 nd state. 

Next we add a non-magnetic impurity scattering center 
to the ellipse at point r , which we model by a local delta- 
function potential 



V{t) = US(t-t ). 



(23) 



The electronic scattering processes are thus described by 
the scattering T-matrix 



r(r„,e) 



U 



l-t/Gg et (r ,r ,e) 



(24) 



The electronic propagator in the presence of the impurity 
potential follows then as 

G(r,r', e ) = G£ e '(r,r',e) 

+ G r et (v, r , e)T(r , e)G r et (r , r', e) , (25) 

and the modified LDOS is 



JV(r,e) = --ImG(r,r,e). 

7T 



(26) 



In our subsequent analysis we will choose selected points 
for the position of the impurity and explore the conse- 
quences for the LDOS in the entire area of the elliptical 
corral. 



IV. RESULTS 

As mentioned above, Fig. 1 shows the free LDOS at 
the energy e = e 42 and for 5 = 0.77E n . The contribution 
of the 43 rd state is very strong along the semiminor axis 
and leads to a sequence of local maxima and minima. On 
the other hand the contributions of the 42 nd state are 
responsible for the structure of the LDOS in the other 
parts of the corral especially near the foci. 

Fig. 2 shows the LDOS when the impurity is placed at 
the right focus; the position of the impurity is indicated 
in the figure by a black bullet. With increasing potential 



strength U one observes that the LDOS is reduced essen- 
tially everywhere except for the positions of the minima 
of the free LDOS and the maxima along the semiminor 
axis, where it remains almost constant. The evolution of 
the LDOS with increasing U is shown in Fig. 5 for a cut 
along the semimajor axis. At the impurity-free focus the 
LDOS continuously decreases with increasing U. Impor- 
tantly, the LDOS is almost symmetrically suppressed at 
the other impurity-free focus; this phenomenon can be 
ascribed to the previously studied so-called quantum mi- 
rage effect. Here we use this acronym to underline that 
the LDOS at both focal points is changing in the same 
manner. We emphasize again that the ridge structure 
with its high maxima along the semiminor axis is con- 
served. So, although the impurity is only at the right 
focus, the semiminor and the semimajor axis still appear 
like symmetry axes - at least within the resolution of the 
chosen color grid. This is the quantum mirage. A similar 
mirage effect persists also, when the impurity is moved 
slightly away from the focus [5] . 
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FIG. 5. LDOS N(r, e) for a cut along the semimajor axis; 
the impurity is at the right focus. The vertical dashed lines 
serve as a guide to the eye to compare the LDOS at the im- 
purity with the LDOS at the mirror position. 

A very different behavior of the LDOS is observed, if 
the impurity is placed at the position of the first maxi- 
mum of the free LDOS on the semimajor axis near the 
center of the ellipse. This is shown in Fig. 3, and the cor- 
responding cut along the semimajor axis is shown in Fig. 
6. Again, as expected, the LDOS in the presence of the 
impurity is reduced with increasing U almost everywhere 
inside the corral. But, surprisingly, the opposite effect is 
observed at the mirror image point of the impurity posi- 
tion with respect to a reflection at the origin. There the 
impurity causes in fact an enhancement of the LDOS. 
Thus, in this case the interference pattern for the per- 
turbed electronic wavefunctions leads just to the opposite 
effect as in the above discussed quantum mirage in Fig. 
2; we therefore call this observation an anti-mirage effect. 
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Note that in the rest of the corral the shape of the im- 
purity LDOS remains very symmetric, but with a slight 
overall enhancement of the LDOS in the left part of the 
ellipse. Similar as in Fig. 2 the maxima of the ridge struc- 
ture along the semiminor axis persist almost unchanged. 
One can observe the anti-mirage effect also when the im- 
purity is put somewhere between ro = (0.16a, 0) and the 
second in i in inn in i it r = (0.33a, 0). 




FIG. 6. LDOS N(r, £42) for a cut along the semimajor axis 
with the impurity at the position r = (0.22a, 0). 

In a third example we place the impurity at the second 
maximum along the semiminor axis away from the center 
(see Fig. 4 and Fig. 7) . This particular impurity position 
has the remarkable consequence that the ridge structure 
along the semiminor axis is entirely wiped out, and what 
is left is very similar to the probability density of the 
42 nd state. This phenomenon occurs for every impurity 
position at any of the local maxima along the semiminor 
axis. For example, if the impurity is placed at the points 
with the highest LDOSs of the elliptic corral, the LDOS 
reacts most sensitively. This in itself may not appear as a 
surprise, because a local perturbation at the space point 
with the highest probability density should indeed lead to 
severe changes in the local electronic structure. Yet, the 
total wipeout of the ridges comes as a spectacular sur- 
prise. If the impurity is placed at a minimum along the 
semiminor axis, the resulting LDOS appears essentially 
unaffected. 

We note that the LDOS strongly depends on the elec- 
tronic density in the corral, i.e. the spatial structure 
of the eigenstates at or close to the Fermi energy. Dif- 
ferent situations can be realized by varying the length 
of the semimajor axis at fixed electronic density. They 
can be physically achieved by studying different corral 
sizes on identical substrates. For example, the energies 
of the 42 nd and the 43 rd state are very close to each 
other so that the LDOSs with the same broadening 8 at 
the eigenenergies €42 and €43 are nearly the same. On the 
other hand, the LDOSs at energies £41 and £44 have a sig- 



nificantly different structure. In particular, the LDOSs 
N(r, £41) and N(r, e 44 ) nearly vanish at the foci. The 
quantum mirage effects as arising from placing impurities 
near the foci, are in this case essentially unobservablc. In 
this sense, the size of the corrals of Co atoms as realized 
in the experiments was a lucky choice for detecting the 
mirage phenomena. 




FIG. 7. LDOS N(r, £42) for a cut along the semiminor axis 
with the impurity at ro = (0, 0.33a). 



V. SUMMARY 

We have identified three distinctly different inter- 
ference phenomena for non-magnetic impurity induced 
changes in the LDOS of non-interacting electrons in el- 
liptic hard-wall quantum corrals. Mirage or anti-mirage 
effects occur depending on the impurity position inside 
the corral. Even an almost complete suppression of pro- 
nounced rich structures in the LDOS can be achieved for 
special choices of the impurity location. These surpris- 
ing phenomena arc a manifestation of quantum mechan- 
ical interference effects which may be tested experimen- 
tally with the already existing elliptic quantum corrals 
on metallic surfaces. 
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